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Cosmological data have provided new constraints on the number of neutrino species and the neu- 
trino mass. However these constraints depend on assumptions related to the underlying cosmology. 
Since a correlation is expected between the number of effective neutrinos N e ff, the neutrino mass 
52 m„, and the curvature of the universe Qk, it is useful to investigate the current constraints in the 
framework of a non-flat universe. In this paper we update the constraints on neutrino parameters 
by making use of the latest cosmic microwave background (CMB) data from the ACT and SPT 
experiments and consider the possibility of a universe with non-zero curvature. We first place new 
constraints on N eff and fi fe , with N eff = 4.03 ± 0.45 and 10 3 Q k = -4.46 ± 5.24. Thus, even when 
Slfc is allowed to vary, N e ff = 3 is still disfavored with 95% confidence. We then investigate the 
correlation between neutrino mass and curvature that shifts the 95% upper limit of ^ m v < 0.45 
eV to ~}2 m v < 0.95 eV. Thus, the impact of assuming flatness in neutrino cosmology is significant 
and an essential consideration with future experiments. 



I. INTRODUCTION 



II. THE EFFECTIVE NEUTRINO NUMBER 



Throughout the previous decades experimental cos- 
mology has benefited from accurate measurements of the 
cosmic microwave background (CMB). The data have 
determined constraints on several cosmological param- 
eters to remarkable accuracy and the ability to constrain 
new physics with the CMB continues to improve. Future 
CMB experiments might even be able to measure B-mode 
polarization and distinguish between neutrino hierarchy 
models. However, when constraining new parameters one 
must be careful when constraints depend on assumptions 
about the underlying cosmology. For example, a correla- 
tion between the neutrino properties and the curvature 
of the universe is clearly expected since, a higher number 
of neutrino species or large mass would introduce pre- 
recombination effects, shifting the positions of the peaks 
in the angular CMB spectrum (cf. [IJ[2]). 

Here we present an update on the constraints of the 
number of neutrino species N e ff and the sum of neu- 
trino masses Sm„ in the framework of non-flat uni- 
verses with f2fc 7^ combining the Wilkinson Microwave 
Anisotropy Probe (WMAP) 7-year [3], South Pole Tele- 
scope (SPT) [3] and Atacama Cosmology Telescope 
(ACT) [5] datasets. 

The paper is organized as follows. In Section|H]we give 
theoretical arguments for why the N e ff and fl k param- 
eters should be correlated. In Section 11111 we discuss our 
method of constraining the parameters N e ff, Em„, and 
fifc . We present the results of the analysis in Section |IV| 
Finally, in Section [V] we conclude and discuss the impli- 
cations of assuming flatness in neutrino cosmology. 



The effective neutrino number N e ff is defined as the 
contribution of neutrinos to the relativistic degrees of 
freedom g*. In a standard physics scenario the particles 
contributing to the total value of g* ~ 10.75 are electrons, 
three neutrinos (and their antiparticles) , and photons. 
Any extra relativistic degrees of freedom can be parame- 
terized in terms of an excess with respect to the standard 
effective neutrino number N e ff = 3 (which more pre- 
cisely is ~ 3.046 after accounting for QED corrections 
and non- instantaneous decoupling of neutrinos) [51 [TJ. 
The neutrino energy density is: 



,r 7/4\ 4/3 
Pu = N ef f - I — j p 1 , 



(1) 



where p 7 is the energy density of photons. A first effect 
of N e ff is related to the primordial helium abundance 
Yp. Changing N e ff affects the freeze-out temperature 
Tfreeze during BBN and therefore the final neutron to 
proton ratio n n /n p [8]. Larger N e ff means earlier freeze- 
out, larger n n /n p , and larger Y p . 

The effect of N e ff on cosmological observables (e.g. 
CMB anisotropy power spectrum and galaxy power spec- 
trum) is emphasized by the epoch of matter-radiation 
equality a eq . In particular, for what concern the CMB, 
an increase in a eq changes the extent of the early Inte- 
grated Sachs- Wolfe effect. The relation between a eq and 
N e ff is given by equating energy densities: 



Prad Pn 



1 + 0.227A, 



eff 



40484 n m h 2 



(2) 
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This shows a linear relationship for a eq (N e ff), which 
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transfers to the baryon to photon ratio at equality [5]: 



3pb 

30496 n b h 2 a 



1 +0.2277V e// Q b 



1.3276 



(3) 



The presence of baryons in the relativistic cosmic fluid 
slows down the sound speed according to the definition, 



1/^3(1 + R), 



(4) 



and so this quantity is also affected at equality by the 
effective neutrino number. This reflects in the size of 
sound horizon at a generic time r [S]: 



/•, : = / dr' c s (t') , 
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y/T+R + y/R + R 



eq 



(5) 



The last equations come from assuming the Universe is 
matter dominated during recombination. As can be seen, 
the sound horizon depends on N e ff through R eq . 




z(10 3 ) 



FIG. 1. A demonstration of how the sound horizon r a changes 
with the effective neutrino number N e ff under the matter 
dominated approximation given by Eq. [5] 

In recent papers (see e.g. [IOTll2j ) it has been found 
that the number of neutrinos is greater than the stan- 
dard model value at more than 2a. The presence of ad- 
ditional neutrinos can be described by a (3+1) or (3+2) 
model with three active neutrinos and one or two ster- 
ile neutrinos [T3J [H] . More exotic solutions may include 



arguments supporting modified dark energy models [151 — 
I17j . In this work we explore possible overlooked parame- 
ter degeneracies that could still favor the standard model 
without the introduction of new physics. 

In [10] the authors provide qualitative arguments for 
how changing the number of allowed neutrinos affects the 
observed values of parameters. One example which we 
use is the relative dependence of distance measurements 
on the Hubble constant. In fact, the sound horizon at 
recombination scales as r s oc 1/H while the distance a 
photon typically diffuses prior to its last scattering goes 
as rd oc l/yfH. This is significant because the response 
of the radiation relative to matter determines the degree 
of damping prior to recombination. In other words, with 
6 S = r s /DA fixed by observation, the angular diameter 
distance, Da, must also decrease as 1/H which is more 
rapid than r d . Thus, the damping increases according to 
9d = "I'd I Da oc \[H [TU|. If these distances vary accord- 
ing to H then they also vary according to any parameter 
correlated with H. In an open universe with nonzero cur- 
vature the effective neutrino number is slightly reduced. 
The theory confirms this because as stated above 9 S is 
constrained by observation which means if N e ff is re- 
duced and f2fc > then r s and Da both increase (see 
Fig. [TJ). However, if the parameter space favors a closed 
universe then there will appear to be a higher number 
of effective neutrinos. This is one of the primary reasons 
for expecting correlation between N e f f and fi^ . 

We conclude this section by stating our purpose further 
constrain the neutrino mass and investigate the effect of 
curvature on this parameter. Finally, we acknowledge 
that there are many papers on the subject of constraining 
neutrino parameters. A non-exhaustive list of additional 
references includes Refs. IT8H3T1. 



III. ANALYSIS METHOD 

In order to fit cosmological models to data we use 
a modified version of the publicly available CosmoMC 
software package [32] ■ This uses a Monte Carlo Markov 
Chain analysis on calculations of the lensed CMB power 
spectrum made with the CAMB package. Our anal- 
ysis combines the following CMB anisotropy datasets: 
WMAP 7-year [3], SPT [4], and ACT [5]. Including 
BAO+i?o simply means we are using the baryon acous- 
tic oscillation (BAO) data of Percival et al. [33] and im- 
pose a prior on the Hubble parameter based on the last 
Hubble Space Telescope observations [34]. We integrate 
spectral data out to l max = 3000. We sample from the 
following parameters: the baryon ftbh 2 , cold dark mat- 
ter Q c h 2 , and dark fl/y energy densities, the scalar spec- 
tral index n s , the optical depth to reionization r, the 
Hubble parameter Hq, and the amplitude of SZ spec- 
trum Asz- We also consider the effective neutrino num- 
ber N e f f , spatial curvature , and the sum of neutrino 
masses ^ m„. 

Finally, we make decisions specific to the high multi- 
pole mode data. We consider purely adiabatic initial con- 
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Parameter 


WMAP7+BAO+Ho 


WMAP7±SPT 


WMAP7±ACT WMAP7±SPT±ACT 


Ha 
n s 

T 

H (km/s/Mpc) 


2.249 ± 0.054 
0.135 ±0.016 
0.721 ± 0.018 
0.979 ± 0.015 
0.086 ±0.014 
75.1 ±3.4 


2.256 ±0.041 
0.130 ±0.0094 

0.722 ±0.015 
0.9808 ± 0.0122 

0.085 ±0.014 
74.0 ±2.0 


2.235 ± 0.047 
0.137 ± 0.012 
0.714 ± 0.018 
0.982 ±0.013 
0.086 ±0.014 
74.6 ± 2.15 


2.258 ± 0.040 
0.129 ±0.0091 

0.722 ±0.015 
0.9803 ±0.0121 

0.086 ±0.014 
73.9 ± 1.92 




4.34 ± 0.88 


3.91 ±0.43 


4.30 ± 0.58 


3.89 ±0.41 



TABLE I. Summary of matching results from WMAP 7-year [3J, SPT [1], and ACT 0. Note that the analyses are modeled 
by the choice to reproduce SPT results, which produces a smaller value for JV e // than expected for ACT data. All datasets 
include BAO and Hq for improved parameter constraints. The quoted errors are given at the 68% confidence levels (CL). 



ditions. When the background data are taken to small 
enough scales the spectra from infrared source emission 
must be taken into account. The IR spectra is domi- 
nated by Poisson power partially from source emission 
clustering at the smallest scales. Thus, a model for such 
effects must be subtracted out from the CMB power spec- 
tra. The resulting adaptations are representative of the 
considerations made during the process of checking the 
code for consistency with established results. It is also 
important to include a Big Bang nucleosynthesis (BBN) 
consistency check during the sampling in order to pro- 
vide analysis consistent with helium abundance measure- 
ments, as proposed in Refs. 35. 3(i . We remark that the 
ACT collaboration did not include the same BBN con- 
sistency condition used by SPT and our analysis. This 
explains why we find a slightly better constraint on N e f f 
thanRef. [5]. 

Additionally, we also constrain the sum of the neu- 
trino masses T,m u . To do this, we use a top hat prior on 
the fractional contribution of neutrinos to the total mass 
density, /„ = i7„/f2 m <E [0,0.5]. Then we extract J^m„ 
from /„ through the standard relation, 

m v = MVt v h 2 eV = §Ah 2 n m f v eV , (6) 

where Vt w = p1/p cr is the neutrino contribution to the 
energy density. 

IV. RESULTS 

Under the flat Universe scenario the constraint im- 
proves to N e ff = 3.89 ± 0.41 at the 68% confidence level 
(see Table III). This result suggests N e ff = 3 is inconsis- 
tent with the data with ~ 95% confidence. 

We then allow the curvature to vary, to determine how 
assuming flatness affects the constraints on N e ff. Fig. [2] 
demonstrates the correlation between fij, and N e ff, 
which agrees with the prediction from Section [n] Inter- 
estingly, the effect of the additional CMB datasets (ACT 
and SPT) increases the correlation between these param- 
eters with respect to WMAP 7-year data alone. This may 
be due in part to the considerable improvement in N e ff 
whereas the uncertainty in the curvature is not notice- 
ably improved by the addition of small scale anisotropy 
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FIG. 2. Correlation between flk and A r e //. The credi- 
ble intervals are given at the 68% and 95% confidence lev- 
els and the markers indicate the locations of the marginal- 
ized values. WMAP+BAO±-Ho is shown in red while 
WMAP±ACT+SPT±BAO±#o is in blue. Note that the ef- 
fect of adding additional datasets is significant. 



measurements. These results suggest that an open uni- 
verse with fewer neutrinos would look similar to a flat 
universe with more neutrinos. We also note that when 
including N e ft as a free parameter in the ACDM±f2fc 
model, the la constraint of Q k = ~0m23±° ™ll found 
in Ref. [3] does not deteriorate significantly for the same 
combination of datasets (i.e. WMAP±BAO±H )- This 
is due to the presence of the BAO data and the Hq prior 
in the analysis, since both probes are sensitive to the ge- 
ometry of the Universe. Therefore, BAO and Hq help to 
break the degeneracy between N e f f and Qk ■ 

Table [n] provides a summary of parameter values for 
runs where Qk and N e ff vary. Here we find N e ff = 
4.03 ± 0.45 and 10 3 Sl fc = -4.46 ± 5.24 at the 68% confi- 
dence level. Therefore, even when Qk is allowed to vary, 
N e ff — 3 is still disfavored with ~ 95% confidence. Note 
that the increased value for Hq is indicative of the known 
correlation between H and f^. We provide an equiva- 
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Parameter 


WMAP7+7V e// +fi fc 


. . .±ACT±SPT 


lOOfW 


2.26 ±0.056 


2.27 ±0.045 


Q c h 2 


0.136 ±0.0169 


0.129 ±0.00915 


n A 


0.721 ±0.0179 


0.723 ±0.0158 


n s 


0.9837 ±0.0157 


0.9863 ±0.0147 


T 


0.0887 ± 0.0148 


0.0894 ± 0.0149 


Ho (km/s/Mpc) 


74.88 ± 3.40 


73.44 ± 2.03 


Neff 


4.61 ±0.96 


4.03 ±0.45 


w 3 n k 


-4.45 ±5.85 


-4.46 ± 5.24 



TABLE II. Summary of constraints while varying fife and 
N e ff. All datasets include BAO and Ho for improved param- 
eter constraints. Errors are at the 68% CL. See Fig. [2] 



Parameter 


WMAP7±BAO±if 


. . .±ACT±SPT 


lOOfW 


2.26 ±0.053 


2.23 ±0.038 


Q c h 2 


0.112 ±0.0036 


0.111 ±0.0029 


n A 


0.719 ±0.0182 


0.726 ±0.0154 


» s 


0.968 ±0.0124 


0.963 ± 0.0092 


T 


0.0897 ±0.015 


0.0873 ±0.014 


Ho (km/s/Mpc) 


69.2 ± 1.6 


69.9 ± 1.37 




< 0.57 eV 


< 0.45 eV 



TABLE III. Summary of the constraint on the sum of the neu- 
trino masses. All datasets include BAO and Ho for improved 
parameter constraints. Errors are at the 68% CL except for 
^m„, which is quoted as a 95% upper limit. 



lent marginalized contour plot of Hq vs. N e ff to empha- 
size the connection (see Fig. [3]) . 

We now turn to the question of how well the datasets 
are able to constrain ^ m, y 



Table III shows the results 
from WMAP in the first column and the result of adding 
the additional datasets in the final column. Although 
the constraint greatly improves the two sigma limit for 
the masses, this is not enough to favor either the stan- 
dard or inverted hierarchy. However, this is not a sur- 
prise because none of the datasets are sensitive enough 
on their own. Forthcoming data from the Planck exper- 
iment and other future experiments will likely improve 
the mass constraint |37j . 

Finally, we investigate the effect of assuming flatness 
while determining an upper bound on ^m„. We in- 
vestigate two models. The first assumes three degener- 
ate massive neutrinos, while the second allows for addi- 
tional relativistic species accounted by AJV e // > 0. We 
define the correlation coefficient pij as the ratio of the 
off-diagonal term of the covariance matrix to the la 
j , so that for two parameters denoted by i and 



errors OiOj 

j we have p^ = 0{jjoiOj. Figure 4] shows that ^2 
and J\- are strongly correlated with a correlation coef- 
ficient of Pn k J2 m v = f° r both models (AN e ff = 
and AN e ff > 0). Furthermore, the degeneracy consider- 
ably increases the uncertainty in the sum of the neutrino 
masses. In fact, with Qj. ^ the 95% upper limit on 




FIG. 3. Correlation between Ho and N e ff. The credi- 
ble intervals are given at the 68% and 95% confidence lev- 
els and the markers indicate the locations of the marginal- 
ized values. WMAP+BAO±i?o is shown in red while 
WMAP±ACT±SPT±BAO±# is in blue. In this case the 
effect of adding additional datasets is also significant. 



Y] rn v more than doubles with respect to the fiat case: 
with m u < 0.95 eV for the model assuming only three 
massive neutrinos and ^2m u < 1.19 eV for AN e ff > 0. 
The strong correlation between curvature and mass is ex- 
pected because massive neutrinos with m v < 0.3 eV are 
still relativistic until recombination so they act as an ad- 
ditional radiative component. As a consequence the pres- 
ence of such massive neutrinos shifts the time of matter- 
radiation equality a eq . Recall the discussion in Section [TT] 
where in this case lower mass neutrinos roughly corre- 
spond to higher N e ff. Neutrinos also leave an imprint 
on the CMB through the early Integrated Sachs- Wolfe 
effect (c.f. Refs. 38, 39 ) which changes the position of 
acoustic peaks. This effect can be compensated for by a 
change in the geometry of the Universe, which weakens 
the constraints on both ^ m v and . See Table IV for 



a summary of cosmological parameters when curvature 
and massive neutrinos are considered. 



V. CONCLUSION 

The resolution of the high effective neutrino number 
in cosmology remains an open question. However, addi- 
tional neutrinos may be due to parameter degeneracy 
or other issues in statistical analysis rather than new 
physics. The focus of this paper has been an argument 
for correlation between the number of effective neutrinos 
N e ff and the curvature of the Universe Qk, which arises 
from the effect of these parameters on distance measure- 
ments. The qualitative argument is confirmed by a sta- 
tistical analysis of CMB anisotropy measurements using 
CosmoMC. 

In this paper we have shown that there is a correlation 
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FIG. 4. Comparison of the correlation between fife and ^ m v 
under the two ANeff models. The model with three mas- 
sive neutrinos is shown in blue while the model with addi- 
tional relativistic species is in red. Intervals are given at 
the 68% and 95% confidence levels and markers indicate 
the locations of the marginalized values. Datasets include 
WMAP7+ACT+SPT+BAO+J7 . The addition of curvature 
allows "^2 ra v to be more than twice the previous constraint. 



Parameter 


AN eff = 


AN R ff > 


100fti,/i 2 


2.24 ± 0.043 


2.26 ± 0.049 


Q. c h 2 


0.118 ±0.0063 


0.134 ±0.0105 


Oa 


0.711 ±0.0216 


0.703 ± 0.0239 


n s 


0.967 ±0.011 


0.982 ±0.015 


T 


0.0864 ± 0.0144 0.0890 ± 0.0145 


H (km/s/Mpc) 


70.6 ± 1.62 


73.1 ±2.03 




7.52 ± 7.74 


3.46 ±8.69 




< 0.95 eV 


< 1.19 eV 


ANeff 





0.995 ± 0.430 



TABLE IV. Summary of the constraint on the sum of the 
neutrino masses when Qk 7^ 0. AN e ff is an additional rel- 
ativistic contribution after considering 3.046 massive neutri- 
nos. Datasets include WMAP7+ACT±SPT±BAO±# . Er- 
rors are at the 68% CL except for ^ m„, which is quoted as 
a 95% upper limit. 



between N e ff and flk that gets stronger when SPT and 
ACT datasets are added to WMAP alone. However, even 
when Slfe is allowed to vary, N e ff — 3 is still disfavored by 
the data with 95% confidence. Although the correlation 
favors a closed universe with fife < 0, if CMB data were 
to favor open models then the neutrino number would 
decrease as predicted. Perhaps the same element of the 
data that favors a closed universe may also be responsible 
for the trend toward a higher N e ff. More importantly, 
we find a strong correlation between curvature and the 
sum of the neutrino masses. 

Future experiments will provide further insight into 
both N e f f and ^2 f^u [40] . Our results are consistent 
with the current understanding of the data available. 
The strongest constraints on these parameters from the 
statistical analysis assuming a flat universe are N e ff = 
3.89 ± 0.41 and Yl m v < 0.45 eV with 95% confidence 
level using WMAP7±ACT±SPT±BAO±# - The con- 
straints are weakened by degeneracy with the curvature 
parameter fife. However, this still represents the contin- 
ued effort toward significant improvements on parameter 
constraints in cosmology. Although the sum of the neu- 
trino masses is significantly improved from the WMAP 
7-year result of m„ < 0.57 eV, the constraint is far 
from being sensitive enough to rule out one of the mass 
hierarchies. Furthermore, we have shown that the mass 
uncertainty more than doubles when Slfc 7^ 0. Based on 
our results and the estimated quality of data for Planck 
and other experiments, it should be possible to determine 
the existence or nonexistence of sterile radiation to much 
greater confidence in the near future. 
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